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ABSTRACT 

We report multi-frequency measurements of scatter broadening times for nine highly dispersed 
pulsars over a wide frequency range (0.6 - 4.9 GHz). We find the scatter broadening times to be 
larger than expected and to scale with frequency with an average power-law index of 3.44 ± 0.13, 
i.e. significantly less than that expected from standard theories. Such possible discrepancies have 
been predicted very recently by Cordes & Lazio. 

Subject headings: ISM: structure - pulsars: general - scattering 



1. Introduction 

Fluctuations in the Galactic electron density 
distribution are responsible for scintillation and 
scattering of radio signals propagating through the 
interstellar medium (ISM). The resulting observ- 
able phenomena have been studied extensively to 
investigate the nature of the density irregularities 
(e.g. Rickett 1977). Pulsars are particularly useful 
as probes of the medium because of their small an- 
gular diameter and spatial distribution with sam- 
ples on many lines of sight through the Galaxy. 
Scattering of pulsar signals through an irregular 
random ISM causes the signal to arrive from dif- 
ferent, multiple ray paths with different geometric 
lengths, so that a pulse which left the source at 
one instant, arrives at the observer over a typical 
time interval, Tgc, commonly the scatter broaden- 
ing time. Further, along different ray paths the 
radiation acquires random phases which cause in- 
terference in the plane of the observer to produce 
diffraction patterns. The pattern decorrelates over 
a characteristic bandwidth Az^d- 

Both the scatter broadening time and decorre- 
lation bandwidth are strongly dependent on fre- 
quency v, i.e. Tsc oc v~°' and At'd cx v°' , and 
are related to each other by 2TiTsc^Vi = C. The 
"constant" G is of the order unity but changes for 
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different geometries and different models for the 
turbulence wavenumber spectrum, P^^ (e.g. Lam- 
bert & Rickett 1999). Knowledge of the spectrum 
would provide valuable insight into the physics 
of the ISM irregularities. A commonly used de- 
scription for the wavenumber spectrum is a power- 
law model with large range between "inner" and 
"outer" scales, and /t~^, (e.g. Rickett 1977), 
i.e. 

Pr^S^l) = 7-^^^^k7?77 exp 



((?2 + nlfn 



q 



(1) 



where q is the magnitude of the three-dimensional 
wavenumber and C^^ is the strength of fluctuation 
for a given line-of-sight. For Kq q Ki, one ob- 
tains a simple power-law model with a spectral 
index /3, i.e. P„„(g) = C^^q-f^ and a = 2/3/(/3- 2) 
(e.g. Lee & Jokipii 1975). For a pure Kolmogorov 
spectrum, /3 = 11/3, we expect a — 4.4. In turn, 
by measuring a, one can infer (3 and hence de- 
tails about the actual wavenumber spectrum. A 
value of /3 = 4 could, for instance, describe a 
medium with abrupt changes in density caused by 
randomly placed, discrete clouds along the line-of- 
sight (Lambert & Rickett 1999). 

In this Letter, we report measurements of a by 
carefully determining Tsc over a wide frequency 
range for a sample of nine highly dispersed pul- 
sars, selected to be observable at 4.8 GHz and lo- 
cated towards the inner Galaxy. Measurements 
of Tsc can be hindered by possible frequency evo- 
lution of the intrinsic pulse shape but we use re- 
sults from previous studies of pulsars unaffected 
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Scatter Broadening for 9 Pulsars 
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Note. — Cols. (2) and (3) give the Galactic longitude and latitude of each pulsar; col. (4) its Dispersion Measure 
and cols. (5)— (9) the scatter broadening with 3cr-crrors in parenthesis at 0.6, 0.9, 1.4, 1.6 and 2.7 GHz. The median 
of normalized X for all fits is 1.37. Col. (10) gives the spectralindex of scatter broadening and col. (11) the spectral 
index of density irregularities for each pulsar with their Icr— errors. Col. (12) gives the ratio k of Tsc at 1 GHz 
scaled by a and Tgc at 1 GHz predicted by the Taylor & Cordes (1993) model. Note that, with the exception of 
PSR B1820— 14, our measured scatter broadening times are significantly higher than those predicted by the Taylor 
& Cordes (1993) model. 
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Fig. 1. — (Top) Observed profile of 

PSR B1815-14 at 1.4 GHz and best-fit of 
the model profile (see section 2). {Bottom) Con- 
tours of near its global minimum (indicated 
by the plus sign) in the c-Tgc plane for the above 
profile. 



by scattering to estimate these pulse shape effects. 
Although one could attempt to measure Av^ in- 
stead, predicted values are well below 50 Hz at 1 
GHz and cannot be measured with the frequency 
resolution of typical filterbanks. This results in 
a surprisingly small number of published spectral 
indices (see Cordes, Weisberg & Boriakoff 1985, 
Johnston, Nicastro & Koribalski 1998). Conse- 
quently, the data presented here are in fact the 
first ever measurements of a for pulsars with high 
dispersion measure (DM). 

2. Observations and Data cinalysis 

Dual circular polarization observations were 
made both with the 100-m radio telescope at 
Effelsberg (1.4, 2.7 and 4.9 GHz) and the 76-m 
Lovell telescope at Jodrell Bank (0.6, 0.9, 1.4 and 
1.6 GHz). Most of the pulse profiles recorded 
at Effelsberg were obtained with a coherently de- 
dispersing backend, the Effelsberg Berkeley Pulsar 
Processor (Backer ct al. 1997). For pulsars listed 
in Table 1, the total available bandwidth ranges 
from 28 to 45 MHz at 1.4 and 2.7 GHz and up to 
128 MHz at 4.85 GHz. At lower frequencies, we 
often obtain higher sensitivity by employing an in- 
coherent hardware de-disperser consisting of 4 x 60 
channels of 0.667 MHz, providing a bandwidth of 
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40 MHz per polarization. Resulting dispersion 
broadening times ^, ^dm, for the pulsars observed 
lie in a range of 0.8 ms < foM < 1-3 ms. De- 
tailed descriptions of the hardware set-up can be 
found in Seiradakis et al. (1995) and Kramer et 
al. (1999). 

At Jodrell Bank, the pulse profiles were ob- 
tained for each polarization with an incoherent 
hardware de-disperser using filterbanks with vary- 
ing number of channels and bandwidtlis. At 
0.6 GHz we used 32 channels of 0.125 MHz band- 
width (2.0 ms < tr,M < 5.4 ms), at 0.9 GHz 32 
channels of 0.25 MHz each (1.1 ms < iDM < 2.0 ms), 
at 1.4 GHz we employed 32 channels of 1 MHz each 
(1.7 ms < tDM < 3.2 ms), while we employed 8 
channels of 5 MHz bandwidth at 1.6 GHz (0.8 ms 
< ioM < 10.0 ms). Details of the system can be 
found in Gould & Lyne (1998). 

Typical observation times were between 30 - 

60 min, depending on observing frequency, band- 
width and flux density of the individual source. 
Incoming signals were folded with the topocentric 
pulse period and two circular polarizations were 
added to produce total power profiles. The signal- 
to-noise ratios of our profiles range from 13 to 67. 

In order to measure the scatter broadening, 
Tsc, we performed least-squares fits of an artifi- 
cial model profile, P^{t), representing the scat- 
tered intrinsic pulse shape, to the observed pro- 
file. For each pulsar, a template P^{t) is con- 
structed as a sum of Gaussian components fitted 
to an unscattered, high frequency profile using 
a procedure developed by Kramer et al. (1994, 
1999). Usually, the high-frequency profile used is 
observed at 4.9 GHz, except for PSR B1820-14 
and PSR B1815— 14 where we use data taken at 
1.6 GHz and 2.7 GHz, respectively, because of the 
low signal-to-noise of their 4.9 GHz profiles. For 
the given DMs, we can expect the scatter broad- 
ening to be negligible at these frequencies (Bhat- 
tacharya et al. 1992), but we discuss the validity 
of this assumption in more detail later. For each 
lower frequency, P^{t) is the convolution of the 
template P'^it) with the impulse response func- 
tion characterizing the scatter broadening, s{t), 
the dispersion smearing across the filterbank chan- 

^tuM = 8.3 DM Au-gw/i''^ M*, where DM is the dispersion 
measure (pc cm^'*), u is the observing frequency (GHz) 
and Ai/BW the filterbanlc channel bandwidth (MHz). 



nel, d{t), and the instrumental impulse response, 

i{t), 

P^I(f) = pT(t)0s(t)0rf(t)0i(t) (2) 

where (8) denotes convolution (Ramachandran, Mi- 
tra & Deshpande 1997). The rise times of the 
receivers and backends are small enough to con- 
sider the effect of i{t) to be ncghgible, while d{t) 
is a rectangular function of width tuM for inco- 
herent de-dispersion and of zero width for coher- 
ent de-dispersion. Williamson (1972) pointed out 
that multi-path scattering due to both thick and 
thin screens leads to an exponential decay of the 
pulse giving s{t) = cxp(— f/rgc)- In general, s{t) 
can assume different functional forms for differ- 
ent geometries reflecting also different rise times of 
scattered pulses (Lambert & Rickett 1999) which 
could lead to uncertainties in the estimation of Tgc- 
However, the frequency dependence of the scatter 
broadening time, i.e. a, can be expected to remain 
unaffected as we discuss later. 

The best fit of the model pulse shape to that 
observed is obtained by minimizing the normalized 

value, 

1 ^ 

^'= (Ar-4W^ ^[i^°-Pf(a,6,c,r.e)]^ (3) 

I /off i—i 

where P° is the observed profile, fTofi the off- 
pulse rms. and N the number of bins in the pro- 
file. During the fit, we adjust four free parameters 
of the model pulse profile P^, i.e. an amplitude 
scale- factor, a, constant offsets in baseline b, 
and phase, c, and the scatter broadening time Tgc- 
We obtain best-fit values and uncertainties for Tsc 
from the x^-contours in the c-Tsc plane. An ex- 
ample for an observed profile, the best-fit model 
profile and x^-contours for the fit is shown in Fig- 
ure 1. 

3. Factors eiffecting scatter-broadening- 
measurements 

Measuring Tsc can be difficult since the intrinsic 
pulse shape often changes slowly with frequency in 
two distinct ways. First, outer profile components 
often have flatter flux density spectra than central 
components (Rankin 1993, Lyne & Manchester 

* usually close to zero 
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Fig. 2. — Scatter broadening times Tgc with their 3c7-errors as a function of observing frequency v for nine 
pulsars. The lines correspond to the linear fit of the form y = ax + K where y = log(Tsc) in ms and 
X = log(z/) in GHz. The dashed lines arc examples of the expected dependence due to a Kolmogorov 
spectrum, i.e. a = 4.4. The bottom right most panel shows a linear function with an average spectral index 
(a) = 3.44. 
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Fig. 3. Spectral index of scatter broadening a 
as a function of DM for our sample of pulsars as 
well as for earlier measurements from Cordes et al. 
(1985) and Johnston et al.(1998). The dotted hue 
a = 4.4 indicates the spectral index predicted by a 
pure Kolmogorov spectrum. The dashed line a = 
4.0 represents the lowest possible index predicted 
within standard theories of the turbulent medium. 



1988 ). Second, pulse widths decrease with in- 
creasing frequency (Sieber, Reinecke & Wielebin- 
ski 1975). While changes in widths are significant 
below ~1 GHz, at higher frequencies the width 
usually saturates to a constant value (Thorsett 
1991, Xilouris et al. 1996, Mitra & Rankin 2001). 
Both effects, if not carefully accounted for, can 
give rise to inaccurate estimation of Tgc- As we 
use high frequency profiles to create our templates, 
stronger outer components could make the tem- 
plate too wide for the lower frequency profiles. 
Due to this effect, wc may tend to underestimate 
Tsc at lower frequencies. Alternatively, at frequen- 
cies around 1 GHz, it is also possible to mistake an 
evolving trailing profile component for parts of a 
scattering tail. This aspect is however not always 
present and hard to predict. In contrast, virtually 
every pulsar's width is increasing towards lower 
frequencies. The model profile would then be too 
small, overestimating Tsc for frequencies well below 
1 GHz. 

In order to quantify the impact of these possi- 
bly competing effects, wc have simulated "typical" 
pulse profiles and their frequency evolution, made 
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them subject to scatter broadening and applied 
our analysis procedure under the additional pres- 
ence of noise. While confirming the above effects, 
the combined profile evolution hardly affects our 
results, i.e. the derived scatter broadening times 
agree very well with the true values. Width evolu- 
tion docs not change the actual values of the scat- 
ter broadening times significantly. The largest im- 
pact is when unnoticed outer components slowly 
evolve at high frequencies, which in our sample; 
is only seen for PSRs B1820-11 and B1849-K00. 
Our simulations show that increased error bars are 
sufficient to account for these effects, so that we 
conservatively quote 3c7 error bars for all scatter 
broadening times in what follows. We addition- 
ally confirm by determinations of a using scatter 
broadening times measured at the two lowest fre- 
quencies only that our results are unaffected by 
possible changes in pulse profiles. 

4. Results and Discussion 

Table 1 summarizes our measurements of Tsc- 
For PSRs Bf750-24, B1758-23, B1805-20, 
B1815-14 and B1820-11, these can be com- 
pared to those obtained by Clifton et al. (1992) 
at a single frequency. We find our more accurate 
scatter times to be in very good agreement with 
their results. Note that all but one of the scat- 
ter broadening times are factors of several larger 
than predicted by the Taylor & Cordes (1993) 
model. However, this model is poorly constrained 
for distant, heavily scattered pulsars. 

Figure 2 displays the measured frequency de- 
pendence of Tsc from which we determine the spec- 
tral indices a listed in Table 1. Earlier multi- 
frequency observations of pulsars with low DMs 
(see Cordes, Weisberg & Boriakoff 1985, Johnston, 
Nicastro & Koribalski 1998) derived spectral in- 
dices that are consistent with a pure Kolmogorov 
spectrum, i.e. a — 4.4 (Figure 3). In contrast, for 
our highly dispersed pulsars we find all individual 
values of a to be significantly lower than a = 4.4 
and usually even lower than a — 4.0 of a "fi = 4- 
model" (Lambert & Rickett 1999). The average 
spectral index is {a) = 3.44 ±0.13 (see also Figure 
2). 

Deriving the spectral index of the wavenum- 

bcr spectrum from equation (1), /3 = 2a/{a — 2), 
we obtain values that are significantly larger than 



those expected from standard theories such as {3 — 
11/3 or /3 = 4 (Table 1). Interestingly, Cordes 
et al. (1985) predicted departures from the Kol- 
mogorov spectrum such as wavcnumbcr enhance- 
ments to be detected in multi-frequency observa- 
tions of high-DM pulsars. Other observational in- 
consistencies with a simple Kolmogorov spectrum 
have been explained by steeper spectra (/3 > 4, 
e.g. Goodman & Narayan 1985), the influence of 
an inner scale (Coles et al. 1987), or physically 
slightly different models (Lambert & Rickett 1999, 
2000). It is, however, interesting to note that re- 
cent results suggest that observations along lines 
of sight of enhanced scattering at low Galactic lat- 
itudes are described by a Kolmogorov spectrum in 
an extended medium (e.g. Stincbring et al. 2000, 
Lambert & Rickett 2000). While the DMs of the 
pulsars in these studies are typically lower, we ob- 
serve an apparently different behavior. 

Lambert & Rickett (1999) point out that the 
scatter broadening function, s(i), depends to some 
extent on the actual geometry and spectral model, 
although our values show that our profiles can 
be described accurately by the model. In the 
worst case, according to their results, our individ- 
ual scatter broadening times could be consistently 
under- or over-estimated by an amount which is 
typically well covered by our 3a error-bars. While 
Tsc computed for a given frequency may change 
slightly, its spectral dependence, a, shows only 
minor variations well within the quoted uncertain- 
ties when testing it by increasing (or decreasing) 
each Tsc by a fixed percentage. We conclude that 
our weak frequency dependence, a < 4, is real 
and cannot be explained by simple application of 
standard theories. 

Anisotropic irregularities may lead to variations 
in the broadening functions for different frequen- 
cies. This causes essentially the phenomena of 
anomalous scattering in the ISM as discussed re- 
cently by Cordes & Lazio (2001). In particular, 
one usually assumes that the transverse extent of 
the scattering screen is arbitrarily large and that 
the strength of the scattering is uniform across 
the screen. However, if this assumption is relaxed 
and/or pulsar distances are large, several anoma- 
lous effects can emerge that can in principle flatten 
the observed frequency dependence of Tsc by re- 
ducing the amount of scattering apparent at lower 
frequencies. In our case, the observed anoma- 
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lous behaviour could be caused by scattering at 
multiple screens with finite extension transverse 
to our line-of-sight. As a consequence, less radia- 
tion reaches the observer at lower frequencies since 
some of the radiation that would be scattered by 
an infinite screen is now lost. For such scatter- 
ing geometries, one can expect to detect several 
breaks in the Tgc power-law spectra (see e.g. Fig- 
ure 3. of Cordes & Lazio 2001). A finer frequency 
sampling and higher signal-to-noisc ratios may be 
able to reveal these spectral features in future ob- 
servations. 

We are grateful to Jim Cordes, Simon Johnston 
and Barney Rickett for helpful discussions. 

Note added in proof. - The equation /? = 
2a/{a — 2) is only valid for /? < 4. As shown 
in Table 1 this results, however, in values of /? 
larger than 4. For turbulence spectra with /3 > 4, 
a relationship /3 = 6 — 8/q; is valid (Cordes et al., 
1986, ApJ, 310, 737; Romani et al., 1986, MN- 
RAS, 220, 19). Applying this to a listed in Table 
1 results however in values (3 < 4. We conclude 
again that our results are not consistent with stan- 
dard theory. We also note that a scatter broaden- 
ing function s{t) = l/\/t cxp(— f/rsc), which is a 
conceivable alternative shape in the limit of one- 
dimensional scatterers that might pertain to situ- 
ations with large scattering (Cordes 2001, private 
communication), does not appropriately describe 
our scattered profiles, in contrast to s{t) used in 
the Letter. 
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